arXiv:1505.00587vl [astro-ph.SR] 4 May 2015 


Draft version May 5, 2015 

Preprint typeset using DTp;X style emulateapj v. 01/23/15 


QUASI-PERIODIC OSCILLATION OL A CORONAL BRIGHT POINT 

Tanmoy Samanta 1 , Dipankar Banerjee 1 , Hui Tian 2 
1 Indian Institute of Astrophysics, Koramangala, Bangalore 560034, India, e-mail: tsamanta@iiap.res.in 
2 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA. e-mail: hui.tian@cfa.harvard.edu 

Draft version May 5, 2015 

Abstract 

Coronal bright points (BPs) are small-scale luminous features seen in the solar corona. Quasi-periodic bright- 
nings are frequently observed in the BPs and are generally linked with underneath magnetic flux changes. We 
study the dynamics of a BP seen in the coronal hole using the Atmospheric Imaging Assembly (AIA) images, 
the Helioseismic and Magnetic Imager (HMI) magnetogram on board the Solar Dynamics Observatory (SDO) 
and spectroscopic data from the newly launched Interface Region Imaging Spectrograph (IRIS). The detailed 
analysis shows that the BP evolves throughout our observing period along with changes in underlying photo- 
spheric magnetic flux and shows periodic brightnings in different EUV and LUV images. With highest possible 
spectral and spatial resolution of IRIS, we attempted to identify the sources of these oscillations. IRIS sit-and- 
stare observation provided a unique opportunity to study the time evolution of one foot point of the BP as the slit 
position crossed it. We noticed enhanced line profile asymmetry, enhanced line width, intensity enhancements 
and large deviation from the average Doppler shift in the line profiles at specific instances which indicate the 
presence of sudden flows along the line of sight direction. We propose that transition region explosive events 
(EEs) originating from small scale reconnections and the reconnection outflows are affecting the line profiles. 
The correlation between all these parameters is consistent with the repetitive reconnection scenario and could 
explain the quasi-periodic nature of the brightening. 
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1. INTRODUCTION 

Coronal bright points (BPs) are bright dynamical fea¬ 
tures seen in quiet-sun and coronal holes. The dy¬ 
namics and evolution was studied in X-rays and EUV 
wavelengths (Vaianaetal. 1973; Golub et al. 1974, 1976a,b; 
Habbal & Withbroe 1981; Zhang et al. 2001; Tianetal. 
2008a; Zhang et al. 2012; Li et al. 2013). They generally live 
for a few hours to few days and have sizes less than 50". BPs 
are believed to be composed of loops connected locally with 
the photospheric bipolar magnetic fields (Golub et al. 1976b; 
Sheeley & Golub 1979). With recent high resolution EUV 
images it is clear that a BP is not a point or simple loop-like 
structure but looks like a miniature active region with mul¬ 
tiple magnetic poles with several connectivities. Moreover, 
depending on the emergence and cancellation of the magnetic 
polarities, the BPs evolve with time and show a lot of dynam¬ 
ics. Theoretical model argues that the interaction between two 
opposite polarities creates a X-point magnetic reconnection 
that locally heats the corona and produces BPs (Priest et al. 
1994; Parnell et al. 1994). The locations of BPs appear to be 
related to the giant convection cells (McIntosh et al. 2014). 
Zhang et al. (2012, 2014) suggested that small bipolar emerg¬ 
ing magnetic loops might reconnect with an overlying large 
loop or open field lines and produce brightenings. They also 
proposed that BPs might consist of two components, one is 
long-lived smooth component due to gentle quasi-separatrix 
layer (QSL) reconnections and other is quasi-periodic impul¬ 
sive component, called as BP flashes. 

Sheeley & Golub (1979) reported that the BPs evolve with 
a 6 minutes time scale. Several observations in X-ray 
and EUV show periodic variation in the intensity of BPs 
over a broad range of periodicity (Sheeley & Golub 1979; 
Nolteetal. 1979; Habbal & Withbroe 1981; Strong etal. 
1992; Kumar etal. 2011; Tianetal. 2008b; Kariyappa et al. 


2011; Chandrashekhar et al. 2013). Some suggested that 
these oscillations are caused by the leakage of acous¬ 
tic waves (p-modes), which propagate along the mag¬ 
netic flux tubes and converts into magnetoacoustic mode at 
higher atmosphere (Bogdan et al. 2003; Kuridze et al. 2008; 
Srivastava & Dwivedi 2010). Others believe that the inten¬ 
sity oscillations are due to repeated magnetic reconnections. 
(Madjarskaet al. 2003; Ugarte-Urra et al. 2004; Doyle etal. 
2006). 

Several studies have been carried out to understand the pe¬ 
riodic nature, but their origin remains inconclusive. Here, we 
study a BP inside a coronal hole (CH) as seen in the AIA 
EUV coronal images (Lemen et al. 2012) and in HMI mag¬ 
netogram (Schou et al. 2012) on the SDO and simultaneously 
with the newly launched IRIS (De Pontieu et al. 2014). Com¬ 
bining imaging and spectroscopic observations, we study the 
dynamical changes within this bright point and its variability. 
We show that the time variability can be explained in terms of 
a repeated magnetic reconnection scenario. 

2. DATA ANALYSIS AND RESULTS 
2.1. Observation and Data Reduction 

Observational data was obtained from IRIS, AIA and HMI 
instruments form 5:14 UT to 6:34 UT on 11 May 2014. 
We used AIA images centered at 335 A, 193 A, 171 A and 
1600 A. AIA and HMI images have 0.6" pixel size and were 
co-aligned. IRIS data was taken in sit-and-stare mode. It was 
pointing towards a coronal hole (centered at 404", —612"). 
Slit-jaw images (SJI) were available only with the 1330 A fil¬ 
ter. We have used IRIS Level 2 processed data which takes 
care of dark current, flat field and geometrical corrections 
etc. The exposure time and cadence of 1330 A SJI and spec¬ 
tra were 8 sec and 9.6 sec respectively. IRIS have a pixel 
size of 0.166". AIA and HMI data was then co-aligned with 
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2.2. Magnetic field evolution 

Figure 2 (a) shows the HMI line of sight (LOS) magne¬ 
togram. Potential field extrapolation was performed by as¬ 
suming a constant a force-free magnetic field (with a = 0) 
(Nakagawa & Raadu 1972; Alissandrakis 1981). The white 
lines connecting different polarities (> ±20 Gauss) repre¬ 
sent the extrapolated field lines. Strong connectivity between 
primary polarity P (positive) and N (negative) can be seen. 
It agrees well with the intensity images as seen in IRIS SJI 
and AIA images. Our ROI lies well within a CH, primarily 
dominated by negative flux over the entire region. Based on 
the main connectivity with the positive flux (P), we have se¬ 
lected a region marked as a white box in Figure 2 (b) to calcu¬ 
late the magnetic flux. Contours with +20 and -20 Gauss are 


Figure 2. (a): HMI LOS magnetogram and extrapolated potential field lines 
(white lines), (b): On the same magnetogram the red and green contours rep¬ 
resent the field above +20 and -20 Gauss respectively. The magnetic fluxes 
(positive and negative) are calculated within the contours (red and green re¬ 
spectively) inside the rectangular white box. (c): The variation of measured 
positive and negative flux over time. 

drawn in red and green respectively and fluxes are calculated 
within these contours. After calculating the fluxes for each 
time frame, smoothed over 3 frames, light curves (LC) are 
shown in Figure 2 (c). There appears to be a good correlation 
between positive and negative flux with a correlation coeffi¬ 
cient (CC) of 0.68. It indicates that both positive and nega¬ 
tive flux amplitudes vary in a similar manner. Movie 1 (avail¬ 
able online) shows small scale magnetic flux emergence and 
cancellation. It is possible that newly emerging small bipolar 
loops are reconnecting with the overlying pre-existing large 


Figure 1. (a): AIA 193 A image showing a coronal hole. The white box represents our Region Of Interest (ROI), covering a bright point in the coronal hole. 
Zoomed view of the ROI as recorded by different AIA channels are shown in (b-e), (f) shows IRIS 1330 A SJI image and (g) shows HMI line of sight (LOS) 
magnetogram. The vertical black line on each image represents the position of the IRIS slit. Intensity within small white box of (b-f) are used to study oscillation 
properties (see Figure 3). The red horizontal tick marked on IRIS 1330 A image (f) is the location where we study the variation of different line parameters of 
the Si IV 1393.76 A line (see Figure 4). A movie of the ROI is available online (Movie 1). 
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IRIS data. IRIS 1330 A and AIA 1600 A were used for co¬ 
alignment. De-rotation was performed on all AIA and HMI 
images to co-align the data cubes. 

Figure 1 (a) shows the coronal hole in AIA 193 A image. 

The magnetogram (Figure 1 (g)) shows that the CH is dom¬ 
inated by negative polarity magnetic field (colorbar can be 
seen in Figure 2). We made a 48 hours movie with AIA 193 A 
images and HMI magnetograms with 1 hour cadence. It re¬ 
veals that the BP appears around 16:00 UT on 10 May 2014 
with emergence of some positive flux and disappears around 
17:00 UT on 12 May 2014 (with a life time ^38 hours) with 
the complete disappearance of the positive flux. The posi¬ 
tive flux cocentration within this BP over the dominated back¬ 
ground negative flux is the probable reason of the existence of 
the BP. Figure 1 shows the BP as seen in various AIA chan- 
nels (b-e), IRIS 1330 A SJI image (f) and HMI LOS magne- 
togram (g). The vertical black line on each image represents 
the position of the IRIS slit. It clearly shows that the IRIS 
slit is crossing one foot point of the bright point, where the 
magnetic field is positive. 
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loop (P-N) as suggested by (Priest etal. 1994; Zhang etal. 
2012, 2014), which could explain the good correlation. In 
section 2.4, we show signatures of reconnection and their re¬ 
lation with flux changes. 

2.3. Imaging observations 

We have focused on the dynamics of one foot point (posi¬ 
tive polarity P). We study the time evolution of this small re¬ 
gion as seen with simultaneous multi-wavelength images cor¬ 
responding to the transition region (TR) and coronal layers. 
We have computed average intensities inside the small white 
box on top of the foot-point as shown in Figure 1 (b-f). The 
size of the white box in AIA was 8x8 pixel ( (4.8") 2 ) and 

IRIS was 15x15 pixel (~ (2.5") 2 ). AIA box size was selected 
slightly larger to accommodate the loop expansion higher up 
and also to reduce movement effects of the loop if any. The 
smoothed light curves (LCs) over 3 time frames are shown in 
the Figure 3. Left panels show (from top to bottom) the LCs of 
IRIS 1330 A, AIA 1600 A, 171 A and 193 A respectively. We 
have performed wavelet analysis (Torrence & Compo 1998) 
on each LC after removing background trend. We use Mor- 
let function, a complex sine wave modulated by a Gaussian, 
for convolution with the time series in the wavelet transform. 
The global wavelet power spectrum are shown in the middle 
panel of Figure 3. A confidence level of 99% is overplotted by 
dotted white line. The confidence level was set by assuming 
white noise (Torrence & Compo 1998). Measured periods are 
printed on the right panels. The global wavelet power plots 
clearly show the presence of periodicities. A dominant period 
around 8 minutes is present in IRIS 1330 A and AIA 1600, 
171 A LCs. Though the global wavelet plots of AIA 171 A 
and 193 A show the presence of strongest peak ~12 minutes 
there is a weaker power around 8 minutes. 

Now we try to explore the source of this oscilla¬ 
tions. Several observational evidence show a positive 
correlation between the EUV and X-ray emission with 
the underlying photospheric magnetic flux (Pres & Phillips 
1999; Handy & Schrijver 2001; Ugarte-Urra et al. 2004; 
Perez-Suarez et al. 2008). Pres & Phillips (1999) have ob¬ 
served that X-ray and EUV emissions are temporally corre¬ 
lated with the photospheric magnetic flux. They also sug¬ 
gested the possibility of several small “network flares” occur¬ 
ring during the lifetime of these BPs. Chandrashekhar et al. 
(2013) have found a good temporal correlation between mag¬ 
netic flux associated with the foot points and the intensity 
brightening and suggested that the possibility of repeated re¬ 
connection scenario. Note that most of these were imaging 
observations. In the following subsection we study the time 
evolution from the IRIS spectroscopic data which provides 
additional information of the possible sources of these oscil¬ 
lations. IRIS sit-and-stare observation provides an ideal op¬ 
portunity to study the time evolution of this foot point as the 
slit position is crossing it. 

2.4. Spectroscopic Analysis 

At first, a single Gaussian fit was performed on the av¬ 
eraged (over all the pixel along the slit and time) profile of 
photospheric Si I 1401.513 A line for absolute calibration of 
wavelength. Now, for our spectroscopic study, we have se¬ 
lected a position (402.40", —574.52") which corresponds to 
one foot point of the bright point loop system. The position is 
marked by a red tick mark on IRIS 1330 A SJI in Figure 1 (f). 
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Figure 3. In each row the left panel shows the variation of intensities of 
different AIA and IRIS channel. Middle panel shows the global wavelet 
power spectrum. The confidence levels are overplotted with dashed lines. 
The right panels show the significant periods as measured from the global 
wavelet spectrum. Intensity variations of AIA images and IRIS 1330 A SJI 
images are calculated from the region inside the white boxes as shown in the 
Figure 1 (b-f). 


An average over 3 pixels along the slit and a running aver¬ 
age of 3 points along dispersion was applied to the spectra 
to improve signal-to-noise ratio. After that, a single Gaus¬ 
sian fit was applied to each IRIS Si IV 1393.76 A line pro¬ 
file to derive line intensity, Doppler shift, and FWHM of the 
line. To compute the asymmetry in the line profile, we per¬ 
formed Red-Blue (RB) asymmetry analysis (De Pontieu et al. 
2009; Tian et al. 2011b; Martinez-Sykora et al. 2011). A sin¬ 
gle Gaussian fit was performed only in the core of the profile 
to find line centroid (similar to that in Tian et al. (2014)). Red 
and blue wings was then subtracted and normalized to peak 
intensity to construct RB asymmetry profile (in percentage). 
Finally, we have constructed the RB asymmetry light curve 
by taking average over 15-40 Km/s velocity range for each 
profile. Positive and negative values represent enhancements 
in the red and blue wings respectively. 

Now, the variation of all the line parameters with time along 
with HMI positive flux and AIA 335 A intensity are shown in 
Figure 4. Different rows (from top to bottom) correspond to 
HMI positive flux, intensity, deviation from average Doppler 
shift (\V — V avg |), FWHM, absolute RB asymmetry (\RB\) 
of Si IV 1393.76 A line and AIA 335 A intensity respectively. 
Now for the correlation studies between different parameters 
and to find out the periodic nature, we have selected a time 
interval of 16-80 minutes. The initial rapid change in flux af¬ 
fects the power analysis, so we have omitted first 16 minutes 
for this analysis. We use trend subtracted smoothed LCs as 
represented by red lines in Figure 4 for easier comparison. 
The wavelet analysis was applied over the smoothed curve to 
investigate the oscillation properties. Global wavelet power 
spectrum along with 99% confidence level are shown in mid¬ 
dle panels. Observed periodicities are then printed on the right 
side. The power analysis on the line parameters show the 
omnipresence of a strong periodicity around 8 and 13 min¬ 
utes. Correlation coefficients (CC) of all the line parameters 
with the absolute RB asymmetry are printed on each panel. 
We correlate different LCs with RB asymmetry LC as the RB 
asymmetry provides a good measurement of the distortion in 
line profile form Gaussian. To investigate it further, we will 
focus on specific instances and will have a closer look at the 
variation of line profiles. We notice sudden changes in the 
line parameters at certain times. In Figure 4, horizontal grids 
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Figure 4. Similar to Figure 3, top row corresponds to HMI positive flux and 
bottom row corresponds to AIA 335 A intensity and other rows correspond to 
different line parameters as calculated from the Si IV 1393.76 A line profiles. 
\RB\ and \ V — V av g \ represents absolute RB asymmetry and the deviation 
from average Doppler shift respectively, where V av g is the average Doppler 
shift. Red lines are the trend subtracted smoothed curves for duration 16 
to 80 minutes. Si IV 1393.76 A line profiles correspond to the location as 
marked in Figure 1 (f) with the red tick mark on the IRIS 1330 A SJI image. 
The correlation coefficient (CC) between all the LCs with the absolute RB is 
printed on the respective panel. 


represent the particular instances where certain changes oc¬ 
cur. The line profiles at Tl, T2, T3, T4, T5, T6, T7 and T8 
time instances are shown in Figure 5 along with averaged (in 
red) and Gaussian fit (in green) profile. The RB asymmetry 
profiles are shown in the lower panel. It clearly shows strong 
asymmetry in the profiles at those particular instances due to 
existence of some secondary emission component. Movie 2 
(available online) shows the evolution of the line profiles. 

The correlation between intensity, FWHM, | V — V avg | and 
\RB\ can be easily explained in terms of sudden flows along 
the line of sight. In active region boundaries, all the line pa¬ 
rameters coherently changes due to quasi-periodic upflows in 
the medium (De Pontieu & McIntosh 2010; Tian et al. 2011a, 
2012). In our study, we observed both red and blueward 
asymmetries time to time (see Figure 5) which is differ¬ 
ent from AR boundary study where predominantly blue- 
ward asymmetry is reported. Our observations can be ex¬ 
plained in terms of both high-speed upflows and down¬ 
flows. Transition region explosive events (EEs) could be a 
possible explanation (Brueckner & Bartoe 1983; Dere et al. 
1989; Innes et al. 1997; Chae et al. 1998; Teriaca et al. 2004; 
Huang et al. 2014). EEs are believed to result from recon¬ 
nection which will affect the line profiles. Reconnection jets 
(regardless of the direction) usually lead to enhancements at 
the line wing (or wings), which would lead to large RB val¬ 
ues, larger line width, larger intensity. The presence of these 
additional flow components usually also leads to large pertur¬ 
bation of the Doppler shift. This sudden changes in the line 
parameters is consistent with the behavior of bursty EEs. We 
find periodic occurrence of EEs at one footpoint of a coronal 
BP. The observed periodicities at one footpoint of the BP at 
the different temperature channels is similar to the observed 
periodic changes in different line parameters which may in¬ 
dicate that recurring EEs are likely producing the oscillatory 
signal seen in the BP. 


Si IV 1393.76 A Line Profiles 
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Figure 5. The panels show line profiles of the Si IV 1393.76 A line at various 
instances as labeled along with their corresponding RB asymmetry profile 
(from 15 to 40 Km/s) in the bottom. Red lines represent average line profile 
over time. Green lines are the Gaussian fits. A movie (Movie 2) of the line 
profiles is available online. 

Chae et al. (1998) reported that EEs happen preferentially 
in regions with weak and mixed polarity magnetic fluxes. 
They also noticed that majority of EEs occur during the can¬ 
cellation of photospheric magnetic field. In our study, we find 
correlation between RB asymmetry and underneath positive 
flux. We use cross correlation technique to find if the two 
LCs are showing similar variation with some time delay (sim¬ 
ilar to Tian et al. (2012)). It shows that \RB\ LC have a time 
delay (TD) of 1.12 minutes from magnetic flux LC. Hence, It 
can be conjectured that the line profiles are strongly affected 
during magnetic flux cancellation phase. It could be that the 
new emerging flux reconnects with pre-existing field which 
cancels the local flux and creates EEs. The periodic behavior 
might be explained by repeated reconnection. One can also 
notice while carefully looking at the variations, that the RB 
asymmetry do not change randomly, it changes slowly with a 
sharp increase around specific instances and then slowly de¬ 
creases. This matches well with a slow reconnection scenario 
(Wang & Shi 1993). Slow magnetic reconnections occurring 
in lower atmosphere could initiate fast reconnection in the TR 
and then decreases slowly. After some relaxation time it re¬ 
peats. 

We also searched for the signature of heating due to recon¬ 
nections. The AIA 335A channel corresponds to ionization 
temperature of about 2.5 x 10 6 K with a wide temperature 
response. In Figure 4, the 335 A light curve shows increase 
of intensity during the time when the line profiles show sud¬ 
den changes. Due to reconnection, magnetic energy releases 
and heats the medium locally at TR. This localized heating 
close to the reconnection point may increase temperature at 
certain pockets in the TR and AIA may see some of these 
emissions. Hence, depending on the heating, higher tempera¬ 
ture emission can enhance during the reconnection time. The 
AIA 335 A light curves show similar property. We wish to 
address this conjecture in our future work while looking at 
coronal spectral lines. 



3. CONCLUSION 
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We study the dynamics of a bright point within a coronal 
hole using combined imaging, spectroscopic and magnetic 
measurements. We focused our analysis on one foot point 
of a bipolar loop structure. Throughout our observation, both 
positive and negative magnetic flux shows correlated varia¬ 
tions. This may suggest that emerging flux interacts with 
the pre-existing overlying fields, which results in reconnec¬ 
tion and cancellation of the flux at the site of bright point. We 
conjecture that the periodic behavior of the positive flux may 
correspond to repeated reconnections which leads to a series 
of continuous periodic brightenings as seen in different EUV 
and FUV lines. We propose that EEs are created due to X- 
point magnetic reconnection and resultant outflows generally 
affect the line profiles. The presence of the secondary compo¬ 
nent emission in the line profiles confirms that. Furthermore, 
we observed enhanced line profile asymmetry, enhanced line 
width, large deviation from the average Doppler shift at spe¬ 
cific instances. The correlation between all these parameters 
is consistent with the scenario of repetitive alteration of the 
line profile by bursty reconnection outflows. We not only ob¬ 
serve similar periodicities at different line parameters of Si IV 
line but also with AIA channels and they are concurrent in 
time and space and hence seems to be related. During EEs 
(see Fig. 4) we see a corresponding change in the magnetic 
field - this is a one-to-one correspondence and certainly indi¬ 
cates the close relationship between the two. 
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